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Abstract. The analysis of extensive numerical data for the percolation proba- 
bilities of incipient spanning clusters in two dimensional percolation at criticality 
are presented. We developed an effective code for the single-scan version of the 
Hoshen-Kopelman algorithm. We measured the probabilities on the square lattice 
forming samples of rectangular strips with widths from 8 to 256 sites and lengths 
up to 3200 sites. At total of more than 10^"" random numbers are generated for the 
sampling procedure. Our data confirm the proposed exact formulaes for the prob- 
ability exponents conjectured recently on the base of 2D conformal field theory. 
Some preliminary results for 3D percolation are also discussed.|J 



1 Introduction 

Percolation is a quite common phenomena in nature. Among known examples 
there are epidemic diseases of garden trees, forest fires, the process of espresso 
preparation, Ising spins at criticahty, etc. ||l|. Ising spins are correlated with 
probability p — 1 — exp(— 2 J/ZcsT) (J is coupling constant, fcs - Boltzmann 
factor and T - temperature) and clusters of parallel spins percolate at the 
critical temperature Tc M. 

Percolation is the simplest example of a critical phenomenon. Scaling re- 
lations for infinite lattices were developed in the same way as for the Ising 
model . Mapping of the percolation model onto a, q = I state Potts model 
[Ej gave a way to predict critical exponents for percolation in two dimen- 
sions using techniques of Conformal Field Theory (CFT) or of Coulomb gas 
representation. 

In the last few years, new insights into the properties of percolation clus- 
ters on infinite lattices at the critical concentration have been developed. 
First of all, Langlands and coauthors found numerically ^ that the proba- 
bility of a cluster spanning the lattice horizontally is an universal function 
of the aspect ratio only, and does not depend on the lattice symmetry or on 
the type of percolation (see also ||^). 

Second, Aizenman proved that there could be more than one incipient 
spanning cluster in two-dimensional critical percolation g . His results were 
in contradiction with a wide-spread belief that the percolation cluster is 



^ Invited Talk at 12th Annual Workshop on Computer Simulation Studies in Con- 
densed Matter Physics, Athens, GA, March 8-12, 1999 
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unique (see, for example ItI). Numerical results for percolation in a strip 
demonstrated that there could be more than one cluster spanning the short 
direction in a strip [B. Later Aizenman's proposal was supported numerically 
for the case of bond percolation M (see also ilQ] ) . 

Finally, Cardy conjectured an exact form for the probability that k clus- 
ters span a large strip or open cylinder of aspect ratio r, in the limit when k is 
large [nil . His result is in a complete agreement with the Aizenman theorem. 

The assumptions of conformal field theory are not rigorously established, 
and it remains important to perform precise numerical tests of the theory. I 
present here an extensive numerical verification of Cardy's proposal. The nu- 
merical data were obtained in collaboration with Sergey Kosyakov p,0,Rl| . 

In Section 2 some known facts on percolation probabilities are shortly 
reviewed. Spanning probability are discussed in Section 3. Most recent find- 
ings on the multiplicity of incipient spanning clusters (ISC), are discussed 
in Section 4, where our previously published results for the ISC probabili- 
ties in bond percolation are compared with the proposed exact result. The 
computational method is described in the Section 5. New results for the fi- 
nite size corrections of spanning probabilities in a strip geometry with up to 
819200 sites are presented in Sections 6 and 7 for open boundaries and cylin- 
der geometry. Preliminary results for the spanning probabilities in a simple 
cubic critical percolation with the samples with up to 1900544 sites are pre- 
sented in Section 8. All averages were done using lO'* strips in all geometries 
investigated. Discussion and acknowledgements finish the paper. 

2 Percolation 

The simplest model of percolation is bond percolation on the square lattice. 
Let us consider a two dimensional (2D) rectangle with an infinite number 
of sites L in both (vertical and horizontal) directions. The probability for 
the bond to be occupied (or placed, or black) is equal to p, < p < 1. 
Consequently, the probability for the bond to be open (or removed, or white) 
is equal to q = 1 — p. li all bonds are closed (that is p = 1) the left side 
is connected with the right one by closed bonds. That is closed, paths from 
the left boundary to the right boundary exist. Clearly, if all bonds are open 
ones {p = 0) there is no way to go by closed bonds through the rectangle. 
It appears that there is some value of the probability pc dividing these two 
possibilities. The value of critical probability pc is equal to 1/2 for bond 



percolation in two dimensions |13 IJ]. This can be seen from the self-duality 
of the 2D square lattice. 

The dual lattice is constructed placing open bonds across the closed bonds 
on the prime lattice, and by drawing closed dual bonds across the open prime 
bonds. An example configuration is plotted in Fig. n^. It is clear from the 
construction that the number of closed prime bonds is equal to the number 
of open dual bonds. If p is the probability that a bond is closed on the prime 
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Fig. 1. Example of a prime lattice and its dual. Closed bonds of the prime lattice 
are denoted by bold full lines, and those of the dual lattice by bold dashed lines. 
Open prime bonds are plotted with thin full lines. Open dual bonds are not shown. 

lattice, then q = 1 —p is the probabihty that the chosen bond be closed on the 
dual lattice. For a given probability p of closed bonds, let us denote by Ph{p) 
the probability that there is a path of connected closed bonds connecting 
left and right sides of the prime lattice. That is Ph{p) is the probability of a 
cluster spanning the lattice horizontally. By P* (q) we will denote probability 
that a cluster of closed bonds spans the dual lattice vertically for a given 
probability q of closed dual bonds. Two relations follows 



Ph{p)+P:{q) = l 

p.ip) + p;:iq) ^ I. 



(1) 



The prime and dual lattices are identical in this construction. Therefore, 
the percolation threshold should be at pc = qc, i-e. Pc — 1/2. A strong proof 



of that statement is the main subject of Kesten's book 1 14 



For the finite system with square size L x L the probability P^ {p\ L) of 
a horizontal crossing is shown on Fig. as a function of p for several L. 
The slope of this curve grows with the lattice size L as L^/'^ , v being the 
correlation length exponent |^ . Is is clearly visible that P^ {p; L) approaches 
a step function as L tends to infinity. Indeed, there are strong mathematical 
results Mj^ql supporting this picture. 

What is the value of the spanning probability exactly at the percolation 
threshold pcl It was believed for a long time that this value was equal to 
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Fig. 2. Spanning probability nh{p,L) as function of bond occupation probability 
p for several values of lattice sizes L = 10, 20, 30, 40, 50 and 200. The larger lattice 
size the larger derivative of spanning probability. 



the value of the percolation threshold |^. Indeed, this is true for the bond 
percolation on the square lattice. From relations (|l|) for the bond percolation 
on the infinite square lattice it is clear that 

P,,{pc) = P:{qc) = P.{Pc) = P^iqc) = 1/2. (2) 

At the same time, the real space renormalization analysis |lq| are based 
on the relation Phijp) — P for the site percolation. In the last case Pc ~ 
0.592746, which is fairly close to 0.5 and reasonable results for the critical 
exponents were computed. However it was recognized later that the value 
of spanning probability Puipc) does not depend on the type of percolation 
as the percolation threshold does. Moreover it turns out that this value is 
invariant and depends only on the aspect ratio of the rectangle [0,009] . 



3 Incipient spanning probability 

All results described in this section as well as in the following sections are 
obtained for critical percolation. The probability p for the bond (or site) to be 
closed is equal to the corresponding critical value Pc throughout the remaining 
text. The value of pc is not universal and depends on the dimensionality, 
lattice type (square, honeycomb, etc.) and percolation type (site, bond, etc.). 
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The spanning probability appears to be universal and invariant under 
shape transformations of the rectangle, also depends on the boundary condi- 
tions P, p7|jl9| ] (see also the mini review by Stauffer [pO[). 

Let us consider a strip with vertical size L^ ~ L and horizontal size Lh = 
rL, where r — L^/Ly is the aspect ratio. Incipient spanning clusters connect 
two segments on the boundary of a macroscopically large strip. Langlands et 
al. H,M found numerically evidence that the probability that a such cluster 
exists is invariant under conformal transformation. Cardy conjectured an 
exact formula for this probability |17| in terms of hypergeometric functions 

where F and 2F1 are the gamma function and hypergeometric function re- 
spectively and T] = {{1 — fc)/(l -I- k))"^, where k defines the aspect ratio r of 
the rectangle r = K(l — k'^)/2K{k'^) as the ratio of two complete elliptic inte- 
grals K{u). The modulus k is associated with the positions (— 1/fc, —1, 1, l/Zc) 
on the real axes mapped under a Schwartz-Christoffel transformation to the 
vertices Xi of rectangle. Thus, the probability Ph is just the probability that 
there are closed paths (Incipient Spanning Clusters) which connect the in- 
terval (a;l,a;2) with the opposite one (x3, x4). This conjecture of Cardy was 
confirmed numerically by Langlands, et al. ||l9[. 

A more practical form of spanning probability was developed by Ziff as 
series expansion in powers of exp(— 27r r) ||2l|] 

Ph{r)^ 2ic(e-'^''/3-|e-^'^'^/3^.^g-i3W3...)^ ^>1 (3) 

= 1 - 2lc (e-^/3'- - |e-W3r _^ ^g-13^/3r^^ )^ ^ < ^ (4) 

where c = r{l)/r{\f « 0.566046680.... Note, the corrections exp(-27r r) 
to the leading behaviour exp(— 27r r/3) decay very fast with increasing r. 

The probability that at least one cluster spans the square with open 
boundaries in both directions is equal to 1/2. This value could be obtained 
from d) or (|) at r = 1. 

In the case of cylinder the spanning probability along the cylinder is 
not known exactly. Simulations by Hovi and Aharony 122] give the value 
Pft,(l) — 0.6366(8) for the aspect ratio r = 1. 

In a forthcoming paper | |31| we compute the spanning probabilities Ph [r) 
for the aspect ratio r in the interval 0.25 < r < 25 for cylindrical geometry as 
well in the interval 0.25 < r < 12.5 in open geometry with the step Ar = 0.25. 
The quality of these results will be discussed in the following sections. 

4 Coexistence of Incipient Spanning Clusters in 2D 

In the following sections we will also compute the probabilities that at least 
two, and even more clusters simultaneously span the lattice from left to right. 
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It was a comiTLon belief until recently that percolation clusters are unique 
on the 2D lattice. Aizennian proposed quite recently in his lecture at the 
Statistical Physics 19 Conference in China p3[ that the number of Incipient 
Spanning Clusters (ISC) in 2D can be larger than one. Later he proved |p| 
that the probability P{k, r; L) that at least k ISCs span horizontally (that is 
along Lfi) the strip {Lfi,L.i,) = {rL,r) with width Ly and horizontal length 
Lh bounded 

where a and a' are different positive constants expected to merge for an 
infinite lattice Q . Indications of the existence of simultaneous clusters in two- 
dimensional critical percolation in the limit of infinite lattices were found in 
computer simulations by Sen |10(] for site percolation on square lattices with 
helical boundary conditions and in a short strips by Hu and Lin [g| using 
their Monte Carlo histogram method. 

We checked numerically the number of spanning clusters in the critical 
bond percolation model on two-dimensional square lattices [g| . We have de- 
termined the numerical values of probabilities 

F(fc,l) = lim P{k,l;L) 

L — 'oc 

for fc = 1,2, and 3 by means of finite-size scaling. The values are given in 
Table Q for the cases of free boundary conditions and periodic boundary con- 
ditions in the vertical direction (i.e. perpendicular to the spanning direction). 



Table 1. Probability to have at least k horizontally spanning clusters in square with 
free boundaries (FBC) and periodic boundaries (PBC). Bond critical percolation. 



P(fc, 1) - FBC 



0.50002(2) 
0.00658(3) 
0.00000148(21) 



P(fc, 1) - PBC 



0.6365(1) 
0.0020(4) 
0.00000014(5) 



Actually, we use fully symmetric lattices, which are self dual even at finite 
sizes. Therefore, relations (g) holds for our finite lattices (for details, see 0) 
and P(l, 1; L) = 1. In fact, the data in Table || support our observation. Finite 
size corrections for the probabilities of at least 2 incipient spanning clusters 
appear to be of the form P(2, 1; L) = P(2, 1) -I- c/L^, with some constant c, 
i.e. proportional to the inverse square lattice size L. We found that the same 
finite size dependence holds also for P(3, 1; L) |^. 

The values of probabilities in Table tell us, that at criticality there 
is percolation in one direction with the probability equal to 1/2. Among 
such samples and in the case of free boundaries, in about 6 samples in 1000 
there are two different clusters percolated horizontally. Finally, 1 sample in 
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Table 2. Probabilities P{k, 1; L) of at least k incipient spanning clusters on critical 
bond square lattices with size {L, L) and free boundaries. Note, that the values of 
P(2, 1; L) are multiplyed by a factor 10^ and those of P(3, 1; L) by a factor lO". 
For each lattice size L the first row is the probability of horizontal spanning and 
the second one is the probability of vertical spanning. 



L 


P(1,1;L) 


P(2,1;L)- 10^ 


P(3, 1;L) •10'' 


8 


0.50005(5) 


7.657(8) 


3.40(15) 




0.50003(4) 


7.660(8) 


3.98(21) 


12 


0.50002(5) 


7.084(9) 


2.57(14) 




0.49995(5) 


7.070(8) 


2.10(13) 


16 


0.50003(7) 


6.855(9) 


1.97(17) 




0.50002(6) 


6.843(8) 


1.79(19) 


20 


0.49990(6) 


6.742(8) 


1.95(14) 




0.50008(5) 


6.745(9) 


1.72913) 


30 


0.49999(4) 


6.650(8) 


1.52(14) 




0.49996(5) 


6.653(7) 


1.52(12) 


32 


0.49999(5) 


6.648(8) 


1.73(12) 




0.50008(7) 


6.642(8) 


1.56(11) 


64 


0.49992(9) 


6.597(9) 


1.33(13) 




0.49999(6) 


6.602(8) 


1.51(14) 



1000000 (one million) contain at least 3 clusters spanning lattice in the given 
direction. The corresponding probabilities for k > 1 are less in the case of 
periodic boundary conditions, that is for cylindrical geometry. 

The data presented in the Table H appear in good agreement with the 
exact form of the probabilities P{k, r) conjectured by Cardy using methods 
of conformal field theory |Ij . He extend the arguments of his early paper [O 
and determined the exact behaviour in the limit of infinite lattices L — > oo 
for the spanning probabilities P{k, r; L) for large aspect ratios r in the case 
of free boundaries 



27r 1 

lim \nPFBc{k,r; L) oc — 7rk{k ~ -) r, 



(6) 



as r -^ cxD for any k. 

Analogously, for periodic boundaries in the vertical direction (spanning 
along a cylinder) he found 

lim \nPpBc{k,r;L) oc -^{k^ - ^) r, (7) 

L— ►CO O 4 

for fc > 2 and ?■ — > cx), and with the different exponent for k ~ I 

lim lnPpBc{^,r;L) (X --— r, (8) 
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in the limit of large aspect ratio r — > oo. 

From the results (plR) one could calculate ratio 

InPFBcjk,!) ^ k{k-\) 
In PpBcik, I) fc2_i ' 

which is I = 0.8 for fc = 2 and | « .857 for fc == 3. Data from the Table g 
give us the values 0.808(10) and 0.857(20), respectively. This close agreement 
with the asymptotic form for r — > oo may be explained by the observation 
[ pr] that the corrections to the leading behaviour are of order exp(— 27rr) 
(compare with expansion (0)). 

The detailed numerical check of Cardy's results (0-p|) is given in the re- 
maining sections. 

5 Computational method. 

Our simulations can be divided into three major parts: 

Sample choosing. There are two ways to generate a sample. The first one, 
is with the number of closed bonds (or sites) fixed, i.e. we choose the 
sample from a canonical ensemble. This method is especially convenient 
for the bond square percolation, when fully self-dual lattices exist even 
for the finite system size |9|. This eliminates not only boundary effects, 
but also 'fluctuations' of the concentration. 

The second way, is with a fixed probability for a given site (bond) to be 
closed, i.e. sampling from a grand canonical ensemble. This method is 
more convenient for the Ziff hull method, for the transfer-matrix method 
and for the Hoshcn-Kopclman algorithm we use in this work. 
It is known that both methods lead to the same results in the 'thermody- 
namic' limit of infinite lattices [Q. It should be noted, that the computed 
averages coincide for both samples, whereas corresponding dispersions are 
quite different, as was found recently by Vasilyev and author p5|| . 
To generate samples, we generate random numbers as result of XOR 
operations (eXclusive OR) applied to the output of two shift registers 
with large length. Usually, we use the pairs of legs {p, q) with (9689, 471) 
and (4423, 1393), though some other combinations from [Q could be used 
as well. This give us an enormously large period of random numbers and a 
lack of acting correlations. Details could be found in papers |26 27| , p8| , p9[ 
and references therein. 

Decomposition to clusters. We use the Hoshen-Kopelman algorithm [M 
in their original "single scan version". That is we don't hold the whole 
lattice in the memory but only the boundary column we start from, the 
current column and the next column. In this way, the algorithm can 
be named as a "dimensionality reduction" Hoshen-Kopelman algorithm. 
This give us the possibility to fit all data together with the instruction 
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code within the 4 MB cache meniory of Alpha workstations. So, we avoid 
the slowest feature of modern computers, i.e. the slowing down of data 
flow from CPU registers to RAM and back, typical for general simula- 
tions. The resulting speed up of simulations is a factor of about 4. 
Check of spanning. We check spanning with the step Ar in aspect ratio 
equal to 1/4 (and even Ar = 1/8 for the largest width L = 256). Namely, 
the sample generation goes from right to left together with the Hoshen- 
Kopelman labeling. After the successful generation of L/ Ar = L/A se- 
quential columns, we compare the cluster labels on the start column with 
the cluster labels of the current column. The resulting information is 
the indicator functions of the aspect ratio Ii{k,r), where k denotes the 
number of independent spanning clusters and i is the sample number. 
The value of this function at a fixed value of variable r is equal to 1 
(there is spanning of k clusters) or (there is no spanning of k clusters). 
This information is added to an output file. After 10^ scans of the strip, 
the file is written onto the disk and a new file is opened. The process 
was repeated 100 times. So, each of 100 files contain the probabilities 

TTh(k,r;L) — ^0^^ J2i=i ^jC^j'^) for the fixed value of strip width L. Er- 
ror bars are computed over 100 such values of probabilities. It should 
be noted, that the probability Trfi{k,r;L) is the probability of exactly k 
clusters, and we have to compute P{k,r]L) — '^'^f^'!Th{j,r; L) (i.e. sum 
from k to oo) to get spanning probabilities Q as defined in Eqs. (||,|1-||). 

This ends the 'measurement' of probabilities. 

The detailed discussion of the computational methods will be published 
elsewhere pT[ . 

The remaining part of the lecture is devoted to the analysis of the resulting 
data P(fc,r;L) g|l|. 

6 Finite size corrections of spanning probabilities. Free 
boundaries. 

In this section we present an analysis of the probability P{k, r; L) that at 
least k clusters span an open strip of width L at least to a distance r. Free 
boundary conditions are used throughout this section. An analysis of the 
probabilities with cylindrical geometry will be given in the next section. 

We generate 10^ different strips of width Ly = L =8, 12, 16, 20, 24, 28, 
32, 48, 64, 128, 256 sites growing up to the length of Lh^rL^ 32000 sites. 
For the site percolation at pc = 0.592746 considered here about 10^^ random 
numbers were generated using 'XOR' combination of the output of two shift 
registers with the large lags. For some cases we repeat simulations with the 
linear congruential rajid48 generator (compare with [0). Comparison did not 
demonstrate any evidence of systematic errors. 

As results of the simulations, we have data for P{k, r; L). Our next task 
is to determine the ISC exponents. For that, we have to take liniL^oo and 
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determine the slope of In P{k,r) as a function of the aspect ratio r. This is 
the first method. Another one, is to determine the slope of hiP(k, r; L) and 
then take the limit of large strip widths L — > cx). 

Theoretically, both methods should give the same result. In practice the 
two methods gave slightly different values which should be comparable within 
error bars. Coincidence of both results could be a good test on the data 
accuracy. Here we will demonstrate how both methods works for the spanning 
probability P(l,r;L) for which we have exact result (Bh and then we will 
apply both methods to the probabilities of multiple ISCs, for which only 
results for very large aspect ratios r are known (^|8|). 

6.1 First method. 



Fig. 3. Coefficient boir) in the fit to the spanning probability in 2D site percolation 
as the function of the aspect ratio r. See, Eqs. ECj. 



First, we have to develop finite size corrections of probabilities. It is 
known from numerics, that the spanning probability for the site percola- 
tion on square lattices with aspect ratio r — 1 behaves at the critical point 
asi 



P(1,1;L)«P(1,1) 



L + Lo 
with some constants 6o = 0.319 and Lq = 1.6. 



(9) 
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Our data are in agreement with this result. We found that ZifF's proposal 
could be generalized | |31| for any aspect ratio r 

P(l,r;£)^P(l,r)+ y;^ (10) 

L + Loir) 

where bQ[r) and Lo{r) are now some functions of the aspect ratio r. LQ{r) 
is a monotonic function vanishing very rapidly and can be set to zero for 
the aspect ratio r ^ 2.5. Function &o(^) is shown on Fig ^. It is interesting 
that the spanning probability for the site percolation does not depend on the 
lattice size for the aspect ratio close to 1/2. We could suggest that this aspect 
ratio r = 1/2 be used for the computation of percolation probability pc with 
a higher accuracy than the usual simulations at r = 1. 

Our resulting function P{l,r) coincides with the exact one (||) within 
error bars. The linear fit to the logarithm of the resulting function lnP(l, r) 
in the interval 1.75 < r < 5.5 gives the slop ei = —1.0476(3) which is very 
close to Cardy's theoretical prediction — vr/S = —1.047197. To get an idea of 
the accuracy, the linear fit to Ziff's approximation in the same interval was 
done and yielded the slope —1.0473(1). So, our first method produces quite 
accurate data. 



6.2 Second method. 

Let us now compute first the slope ei{L) of the curves lnP(l,r;i) with 
respect to the aspect ratio r for each strip width L, i.e. ei{L) and then take 
the limit 



lim ei(i) -^ ei. 

L — ^oo 

Values of ei(i) are given in the Table 0. A fit to the data in the table give 

with the value of ei equal to the exact one in 4 digits. 

It is interesting that the fit of values ei(L) with some exponent 6, i.e. with 
ei(L) = ei+A/L« yield the result ei(L) = -1.05093(1) + 0.8365(2)/LO-8522(i) 
with too large value of ei and the exponent 9 — 0.8522 close to the irrelevant 
exponent, introduced by Stauffer (see, f.e. |^) and included in the correc- 
tions to the scaling in the detailed analysis of invariance in two dimensional 
percolation [^. Our data does not support this exponent but rather the 
expansion in powers of 1/L which could be effectively written as the term 
A/{L + Lq). Such a term was used first in a paper by Levinstein, et al. (see, 
f.e. p31)n in their first-time computation of the correlation length exponent 

^ We are glad to acknowledge discussion of this point with E.I. Rashba and 
A.L. Efros, who stressed that they was introduced this term "to make the line 
as straight as possible". 
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and successfully used later in the context of spanning probabilities by Ziff 
in his computation of the percolation threshold value Pc = 0.592746 for the 
site percolation, which is the most precise value (see also p3). We found this 
form to be very helpful in the data analysis. 



Table 3. Slope ei{L) of the lnP(l,r; L) as a function of the aspect ratio r = Lh/Ly 
for several strip widths L„ = L. Slopes are computed in the range 1.75 <r< 5.5. 



L 


fiW 


a 


8 


-0.909926 


.0001 


12 


-0.949355 


.00003 


16 


-0.971482 


.0001 


20 


-0.985734 


.0001 


24 


-0.995344 


.00005 


28 


-1.001989 


.00008 


32 


-1.007151 


.0001 


48 


-1.019926 


.00004 


64 


-1.026430 


.00007 


128 


-1.036904 


.00007 



oo -1.047032 .00002 



-0.2 



-0.4 



-0.6 



-I 



Ai- fi'm 



Fig. 4. Derivative of the logarithm of the spanning probability P(l, r; L) as a func- 
tion of the aspect ratio r for the strip width L = 8, 12, 16, 24, 32, 48, 64, 128, 256 
from top to bottom. Site percolation on a square open strip. 
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The second method could be modified in the following way. Taking the 
derivative of hiP(k, r) with respect to aspect ratio r we get 



— - lnP(fc, r) oc 
dr 



3 ^ 2' 



e/c- 



(12) 



This give us the idea the to compute e^ (r, L) taking the numerical deriva- 
tive of our data In P{k, r; L). The result is shown on Fig. ^ for site percolation 
in a strip with several widths L. The slope of lnP(l, r; L) is constant in some 
region of r. Deviation below the left-hand value Tmin ~ 1-5 is due to the cor- 
rections exp(— 27rr) to the leading behavior as seen from Ziff 's approximation 
(0) and becomes negligibly small starting from r = 1.75. Large fluctuations 
for aspect ratios larger than Tmax ~ 5.5 is due to the small values of proba- 
bilities P(l, r; L) for large r and, therefore, error bars become of the size of 
the fluctuations (not shown on the Fig. H for the reason of picture clarity). 
The Fig. H shows a magnified picture with the same derivative for the values 
of the aspect ratio 1.75 < r < 5.5 together with the error bars. 
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Fig. 5. The same as Fig.^ for the aspect ratio in the region 1.75 < r < 5.5. 
Derivative of the logarithm of the spanning probability P{l,r;L) for strip widths 
L = 8, 12, 16, 24, 32, 48, 64, 128, 256 from top to bottom. 



We now take the simple average of ei(r, L) over r in the mentioned region 
and get values of £i{L). In the same manner as before, we could now take 
ei = limi^oo ei(i) ~ 1.0473(1) which is consistent with the values obtained 
by other methods and coincides well with the exact value. 
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It should be noted that the accuracy of this method is not as good as 
that of the previous ones because we take simple numerical derivatives, which 
magnify all fluctuations. Nevertheless, the method of derivatives give us clear 
idea on the interval (rmimfmax) of aspect ratio over which we can safely 
compute exponents. So, it is helpful as a supplement to the two methods 
described above. 



6.3 Two and three ISC's in an open strip. 

The same analysis of the data for probability of at least two ISC leads to 
the exponent 62 = —6.294(6) which is within two sigma from the exact value 
—2tt = —6.2832. The corresponding exponent for the three ISC's are also 
computed 63 = —15.51(1). This value is well coincides with the exact value 
— 28/37r « —15.708. We got relatively large deviation because the value of 
such probability decays very fast with the aspect ratio r and practically we 
take averages only over the rather small interval of aspect ratio 0.5 < r < 
1.25. The good agreement could be explained by the small corrections of 
order of exp(— 27rr) to the leading behaviour [nil. 



7 Periodic boundary conditions. 

Periodic boundary conditions in vertical directions lead to different values of 
spanning probabilities. The geometry is equivalent to an open ended cylinder 
with the circle length Ly = L and "horizontal" length Lh = tL. Clusters 
can span from the one cylinder end to another one also rotating around the 
cylinder. This leads to larger values of the one cluster spanning probability 
than in the case of open strip described in above. Actually, we obtain a 
spanning probabihty equal to 0.6365(1) which is larger than 0.5 in the case of 
the open strip and in a good comparison with the result of the paper |g^ and 
more accurate. So, the spanning probability depends not only on the aspect 
ratio but on the boundary conditions as well [^jl9,2^. As a consequence, the 



values of multiple ISC probabilities should be smaller than the corresponding 
ones for the open strip. 

We apphed both methods described in the previous section to the data 
analysis of the spanning probability P(l, r; L) along an open ended cylinder. 

Using the first method we compute the slope ei of the logarithm of the 
spanning probability and get ei = —0.6547(2). Fig. || shows the probability 
together with the best fit. 

The second method based on the slope ei(L) calculation for each given 
lattice size L = 8, 12, 16, 20, 24, 28, 32, 36, 40, 48, 64 and then taking the limit 
of infinite lattice, finally leads to the value ei = —0.65448(5). Both values 
are very close to the exact value — 57r/24 = —0.654498. The second method 
produce much better accuracy. 
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Fig. 6. Spanning probability P(l,r) along cylinder as function of aspect ratio r 
denoted by open rhombs with error bars. Solid line is the best fit to probability in 
the interval 1.75 < r < 5.5. 



Finite size corrections to the limiting values of the slope are similar to 
those for free boundaries. The only difference is in the finite size corrections 
of ei(i). The leading term, proportional to the inverse width vanishes, and 
the second power plays the main role, namely ei(L) = ei — 0.309(6)/L^. (We 
recall, that in the case of open boundaries the finite size corrections were 

ei(L) = ei + c/(L + Lo).) 

The l/L"^ finite size corrections also hold for e2{L) and e3(i), whose values 
are e2{L) = —7.852(1) and e^[L) = —18.11(15), also quite close to the exact 
values -57r/2 = -7.8540 and -357r/6 = -18.33. 

8 Incipient spanning clusters in simple cubic lattice 

Aizenman, in his seminal paper ||^, also proposed the ISC exponents in di- 
mension d from lower d; = 2 to upper d„ = 6 critical dimensions 

Ph{k,r) oc exp(— afc^i^^) = exp(— e^r). (13) 

The preliminary analysis of the data obtained for the probabilities of 
clusters spanning between two opposite planes of the open cube (rL, L, L) 
at the critical site percolation threshold pc = 0.3116 [ p4|j35| ] are presented 
in the Table ^. We compute values of ek using the second method described 
above for system sizes L = 8, 16, 32, 48, 64 and aspect ratio up to r = 7.25, so 
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some samples contain 1900544 sites. The averages were taken over 10^ such 
samples. 

A possible fit to the data in Table H is a/c ex -7r(fc^-^-0.22fc + 0. 11), which 
supports Aizenman proposal. Wc cannot exclude, nevertheless, the power in 
exponent being less than d/{d — 1) = 1.5. 

In their preprint |36[| , Sen and Bhattacharjee proposed on the basis of their 
simulations that in the directed 3D percolation the ISC exponent depends 
with the same power 2 on the number of clusters k. Here we compute directly 
exponents as the slope of the logarithm of the corresponding probabilities, 
whereas in the Sen and Bhattacharjee preprint only values of probabilities for 
one fixed aspect ratio r — 1 were determined. We are not sure their analysis 
could distinguish the powers 1.5 and 2. We could only say that analysing our 
very precise data, we believe the value of power to be not larger than 1.5. 

More detailed simulations and analysis are in progress. 



Table 4. Exponents et depending on the number k of ISC clusters in simple cubic 



site critical percolation. 



k 


Efc 


a 


1 


-2.7641 


0.0004 


2 


-7.935 


0.001 


3 


-14.68 


0.01 


4 


-21.7 


0.2 



9 Discussion 



We compute the values of the new exponents connected with the probabilities 
of a number of incipient spanning clusters in two dimensional percolation at 
criticality. 

The results of comparison of the computed values with the exact ones 
are summarized in Table ^. Our data strongly support the exact results (|§-||) 
for the number of incipient spanning clusters in two dimensional percola- 
tion developed by Cardy in [f^Ol ^'^'^' therefore, the validity of conformal 
invariance for multiple spanning probabilities. 

Finite size corrections of the data computed with a high accuracy demon- 
strate that at criticality, corrections to all quantities analyzed depends on 
the inverse system size or on the inverse square size. We did not find any 
evidence for the irrelevant exponent 9i = 0.85 as discussed in p2[ . 

We stress the possibility of self duality for finite square lattices in the case 
of bond percolation. 

We generalize finite size corrections developed by Ziff for the site perco- 
lation in a square geometry to any aspect ratio. We found that finite size 
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Table 5. Comparison of computed data with the exact results. Site percolation in 
two dimensional square lattice with free boundaries (FBC) and periodic boundaries 
(PBC). 



k 


Free boundaries 


Periodic boundaries 


Ok Exact Ok Comp. 


Qfe Exact Ofe Comp. 


1 

2 
3 


-f = -1.047197 -1.0473(2) 
-2-K = -6.28318 -6.294(6) 
-57r = -15.70796... -15.51(1) 


-|f = -0.654498 -0.65448(5) 
-^ = -7.8540 -7.852(1) 
_352r ^ _ig 33 -18.11(15) 



scaling of the probability of cluster spanning along a strip with open bound- 
aries is described by a nonmonotonic function ba{r) which vanishes at the 
aspect ratio r ~ 1/2, or nearby. 

Our prehniinary data support Aizennian's proposal for the ISC exponent 
in three dimensional critical percolation. More data and analysis are necessary 
for definite check. 
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